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Abstract. In this talk I discuss the role of proto-globular clusters as the dom- 
inant sources of radiation that reionised hydrogen in the intergalactic medium 
(IGM) at redshift z ^ Observations at lower redshift indicate that only a 
small fraction, (fesc), of hydrogen ionising radiation emitted from massive stars 
can escape unabsorbed by the galaxy into the IGM. High redshift galaxies are 
expected to be more compact and gas rich than present day galaxies, conse- 
ly^ ' quently {fesc) from their disks or spheroids might have been very small. But if 

the sites of star formation in the galaxies are off-centre and if the star formation 
efficiency of the proto-clusters is high, then the mean {fesc) calculated for these 
objects only, is expected to be close to unity. 
\ Here I argue that this mode of star formation is consistent with several 

^n , models for globular clusters formation. Using simple arguments based on the 

observed number of globular cluster systems in the local universe and assuming 
that the oldest globular clusters formed before reionisation and had {fesc) '^1,1 
show that they produced enough ionising photons to reionise the IGM at z ^ 6. 
I also emphasise that globular cluster formation might have been the dominant 
mode of star formation at redshifts from 6 to 12. 



1. Introduction 



C/3 . 

. In this talk, using simple arguments, I emphasise the important cosmological role 

of globular cluster (GC) formation at high-redshift. I show that the formation 
of GCs may have been the dominant mode of star formation near the epoch 
^ ' of reionisation and have contributed signific antly to it. Th e material presented 

in this talk is based on published work by iRicotti Sz Shull (,2000. ) and .Ricottil 
(|i002). 

^Observation of Lya absorption systems toward high-redshift quasars (Beck er et alJ 

I2OOII ) indicate that the redshift of reionisation of the interg alactic medium (IG M) 
is Zrei ~ 6. The recent result from the WMAP satellite (iKogut et al.ll20n,tl of 



an earl y epoch of reionisation will not be addressed in this talk. The reader can 
refer to IRicotti Qstrikeil (|2003al lbl) if interested in this topic. 

A key ingredient in determining the effectiveness by which galaxies pho- 
toionise the surrounding IGM is the parameter {fesc), defined here as the mean 
fraction of ionising photons escaping from galaxy halos into the IGM. Cosmo- 
logical simulations and semi-analytical models of IGM reionisation by stellar 
sources find that, in order to reionise the IGM by 2: = 6 — 7, the escape frac- 
tion from galaxies must be relatively large: {fesc) ~ 10% assuming a Salpeter 
initial mass function (IMF) and the standard ACDM cosmological model. The 
assumption of a universal star formation efficiency (SFE) is consistent with the 
observed values of the star formation rate (SFR) at < z < 5 and total star 
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budget at z = 0. However, the assumption of a constant (fesc) does not seem 
to be consistent with observations. An escape fraction {fesc) ~ 1 is required 
for reionisation at z ^ 6 but the ioni sing background at z ~ 3 is consistent 
with {fesc) ~ 10% (|Bianchi et al.ll2nni] ). Small values of {fesc) at 3 are also 
supported by direct o bservations o f the Lyc emission from Lyman-break and 
starburst galaxi es. iGiallongo et alJ (|2002 l find an upper limit {fesc) < 16% at 

3 (but see ISteidel et al.ll2nnl^ and observations of low-redshift starbursts 
are consistent with { fesr) upper limits rang ing from a few percent up to 10% 
( Hurwitz et J]ll997l : fPeharveng et al.ll2nni h. 

Theoretical models fe-o.- IPov^taL II2OOO) for the radiative transfer of ion- 
ising radiation through the disk layer of spiral galaxies similar to the Milky 
Way find {fesc) ~ 6%. At high redshift the mean value of {fesr) is expected 
to decrease almost e xponentially with increasing redshift (lEicotti fc Shullll2nnnl : 
IWood fc Loeio ''20n0'1: at z > 6, {fesc) ;S 0.1 - 1% even assuming star formation 
rates typical of starburst galaxies. The majority of photons that escape the halo 
come from the most luminous OB a s sociat ions located in the outermost parts of 
the galaxy. Indeed iRicotti &: Shulll ()200(]| ) have shown that changing the lumi- 
nosity function of the OB association and the density distribution of the stars 
has major effects on {fesc) (see their Figs. 8 and 9). 

A star formation mode, in which very luminous OB associations form in the 
outer parts of galaxy halos, may explain the large {fesc) required for reionisation. 
Globular clusters are possible observable relics of such a star formation mode. 
Their age is compatible with formation at reionisation or earlier. Because of 
their large star density they survived tidal destruction and represent the most 
luminous tail of the luminosity distribution of old OB associations. In § 123 I 
show that several models for the formation of proto-GCs imply an {fesc) ~ 1. I 
will also show that the total amount of stars in GCs observed today is sufficient 
to reionise the universe at z ~ 6 if their {fesc) ~ 1- This conclusion is reinforced 
if the GCs we observe today are only a fraction, 1/ fdi, of primordial GCs as a 
consequence of mass segregation and tidal stripping. 

Ill § ED I briefly review a few observational properties and in § ED formation 
theories of GCs that motivate the assumption of {fesc) ^ 1; in §ElI discuss the 
model assumptions in light of GC observations and present the results. In § El 
I present my conclusions. 



2. Condensed review on GC systems 

Most galaxies have a bimodal GCs distribution indicating that luminous galaxies 
experience at least two major episodes of GCs formation. The bulk of the 
globulars in the main body of the Galactic halo appear to have formed during a 
short-lived burst (~ 0.5 — 2 Gyr) that took place about 13 Gyrs ago. This was 
followed by a second burst associated with the formation of the galactic bulges. 

The method for determining the absolute age of GCs is based on fitting 
the observed colour- magnitude diagram with theoretical evolutionary tracks. 
The systematics in the evolutionary model and the determination of the cluster 
distance are the major sources of errors. Rec ent determinations of the absolute 
age of old GCs find tgc = 12.5 ± 1.2 Gyr ()Chabover et al.l[T99§ ). consistent 



with radioactive dating of a very metal-poor star in the halo of our galaxy 
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( Cavrel et all 120011 '!. Relative ages of Galactic GCs can be determined with 



greater accuracy, since many systematic errors can be eliminated. In our Galaxy 
Atgc = 0.5 Gyr, but dif ferences in age betw een GC systems in different galaxies 
could be Atgc ~ 2 Gyr (IStetson et al.lll99fir ). 

The GC specific frequency is defined as the number, N, of GCs per My = 
—15 of parent galaxy light, Sn = N x iqOA{Mv+15) ^jjarris & van den Bergh 
1981). The most striking characteristic is that /Sat (Ellipticals) > S'at (Spirals). 
Sn = 0.5 in Sc/Ir galaxies ( Harris, 1991), S]v = 1 in spirals of types Sa/Sb , 
57V = 2.5 in field ellipticals (iKundu fc Whitmord l200lh . Converting to lumi- 
nosity {Lv/Lq = iO"0-4(^V-4.83)) jjg^^g ^ {L/Lq)Sn/8.55 x 10^. I can 
therefore calculate the efficiency of GC formation defined as. 



Mgc _ fdiNrUgc _ SNfdi 



M M (M/L) 



X 0.00585, (1) 



V 

where Mgc is the total mass of the GC system, rrigc = 5 x 10^ Mq is the mean 
mass of GCs today, M is the stellar mass and {M/L)v is the mass to light 
ratio of the galaxy. In the next paragraph we show that, because of dynamical 
evolution, mgc and N are expected to be larger at the time of GC formation than 
today. Therefore, the parameter fdi > 1 is introduced to account for dynamical 
disruption of GCs during their lifetime. 

The IMF of GCs is not known. The present mass function is known only 
between 0.2 and 0.8 Mq since high-mass stars are lost because of two-body 
relaxation and stellar evolution processes. Theoretical models show that the 
shape is consistent with a Salpeter-like IMF. The mean metallicity of old GCs is 
Z ~ 0.03 Zq. One of the most remarkable properties of GCs is the uniformity 
of their internal metallicity A.[Fe/H] ^0.1. This implies that the bulk of the 
stars that constitute a GC formed in a single monolithic burst of star formation. 
A typical GC emits S" ~ 3 x 10^^ s~^ ionising photons in a burst lasting 4 Myrs: 
about 300 times the ionising luminosity of largest OB associations in our Galaxy. 

During their lifetime GCs lose a large part of their initial mass or are com- 
pletely destroyed by internal and external processes. Numerical simulations 
show that about 50%-90% of the mass of GCs is lost due to external processes, 
depending on the host galaxy environment, initial concentration a nd IMF of the 
proto-GCs (|Chernoff fc WeinberdlT99nl : iGnedin fc OstrikeI^ll997^ . Many of the 
low-metallicity halo field stars in the Milky- Way could be debris of disrupted 
GCs. The mass in stars in the halo is about 100 times the mass in GCs. There- 
fore the parameter /^j, defined in §|T1 could be as large as fdi = 100. Overall fcti 
is not well constrained since it depends on unknown properties of the proto-GCs. 
According to results of N-body simulations fdi should be in the fdi^'^ — 10 
range. 



3. Why is {fesc) ~ 1 plausible for GCs? 

I discuss separately two issues: (i) the (fesc) from the gas cloud in which the 
GC forms, and (ii) the (fesc) through any surrounding gas in the galaxy. 

(i) The evidence for (fesc) ^ 1 comes from the observed properties of present- 
day GCs. The fact that they are compact self-gravitating systems with 
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low and uniform metallicity points to a high efficiency of conversion of gas 
into stars. A longer timescale of star formation would have enriched the 
gas of metals and the mechanical feedback from SN explosions would have 
stopped further star formation. If w 10% of the gas is converted into 
stars in a single burst (with duration < 4 Myr) at t he centre of a s pheroidal 
galaxy, following the simple calculations shown in lRicotti Sz ShuU (200^) 
[see their eq. (18)] at z = 6 we have, fesc = 1 — 0.06(1 — f*)'^/ f* ~ 50%. 

(ii) The justification for {fesc) ~ 1 is model dependent but in general there are 
two main arguments: a) the high efficiency of star formation /^,, and b) 
the sites of proto-GC formation in the outermost parts of the galaxy halo. 

In the "cosmological objects model" (30 < z/ < 7) of iPeebles Dick3 
(1961) GCs form with efficiency /=„ 100%, implying {fesc) = 1 (note 
that such a high is n ot found in numerical simulations of first object 
formation (iRicotti et al.i ■2002a.bi')l. I n "hierarchical formation m odels" 
(10 < Zf < 31 (|Harris fc PudritdllOOi : iMcLaughhn fc Pudritzlll99fil . e.g.,) 
GCs form in the disk or spheroid of galaxies with gas mass Mg ~ 10^ - 10^ 
Mq. Compact GCs survive the accretion by larger galaxies while the rest 
of the galaxy is tidally stripped. Assuming that 1 — 10 GCs form in a galaxy 
with Mg ~ 10^ - 10^ M0 imphes f^ ~ 10% and therefore {fesc) > 50%. 
{fesc) is larger if proto-GCs are located off-centre {e.g., if they form from 
cloud-cloud collisions during the galaxy assembly) or if part of the gas in 
the halo is collisionally ionised as a consequence of the virializ ation process. 
In mo dels such as the "super-shell frapmentation " (zf < 10 ) of Tanigu chi et aP 
( 19991 ) or the "thermal instability" {zj < 7) of lFall fc Re es (1985), {fesc) 
« 1 since proto-GCs form in the outermost part of an already collisionally 
ionised halo. 

In summary, since {fesc) depends strongly on the luminosity of the OB as- 
sociations and on their location, proto-GCs, being several hundred times more 
luminous than Galactic OB associations, should have a comparably larger {fesc)- 

4. Method and Results 

In this section I estimate the number of ionising photons emitted per baryon per 
Hubble time, Mph, by GC formation. In § 14.1.1 1 derive Mph assuming that all 
GCs observed at z = formed in a time period Atgc with constant formation 
rate. In § 14.2.1 1 use the Press-Schechter formalism to model more realistically 
the formation rate of old GCs. 

Table 1. Star census at z=0. 



Type 


(%) 


Sn 


(M/L)v 


(%) 


Sph 


6.5+^ 


2.4 ±0.4 


5.4 ±0.3 


0.26 ±0.06 


Disk 


9+1.5 
^-0.5 


libO.l 


1.82 ±0.4 


0.32 ±0.1 


Irr 

Total 


n ir:+0.15 
'-'■-'^^-0.05 
q+5.5 


0.5 


1.33 ±0.25 


0.22 ±0.04 
0.3 ±0.07 
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4.1. The simplest estimate 

I start by estimating the fraction, ujgc, of cosmic baryons converted into GC stars. 
By definition iVgc = uj^€gc, where uj^, is the fraction of baryons in stars at z = 
and €gc is the efficiency of GC formation defined in § [TJ In ah the calculations I 
assume = 0-04. In T able [H I summarise the star census at z = according to 
iPersic fc Sahi"^ (Il99^ and I d e rive e nr using eq. (QJ, assuming /^j = 1. Using 



similar arguments iMcLaughlinI (|l999l l finds a universal efficiency of globular 
cluster formation Cgc = (0.26 it 0.05)%, in agreement with the simpler estimate 
presented here. It follows that iVgc = fdiC^-^-fj x 10""^) at z = 0. 

The total number of ionising photons per unit time emitted by GCs is 
r]u>gc/ /^tgc, where rj is the number of ionising photons emitted per baryon con- 
verted into stars, and ujgc ~ 2.1ujgc takes into account the mass loss due to stellar 
winds and SN explosions adopting an instantaneous-burst star formation law. 
GCs did not recycle this lost mass since they formed in a single burst of star 
formation, r] depends on the IMF and on the metallicity of the star. I calculate 
r] using a Sal peter IMF and star in etallicitv Z = 0.03 Zq (see § EiJ) with Star- 
burst99 code ( Leitherer et al.lll999l ). and find t] = 8967. The number of ionising 



photons per baryon emitted in a Hubble time at z = 6 is, 

where I have assumed (fesc) = 1 and Hubble time at z = 6 t// = 1 it 0.1 Gyr. I 
expect I < fdi ^ 100 and 0.5 ^ At,, ~ 2 Gyr. A conservative estimate of fdi ~ 2 
and Atgc ^ 2 Gyr (i.e., W < Zf < 3) implies fdi/ Atgc ^ 1. The IGM is reionised 
when Mph = Chii, where Cuu = {''t-'hii) / {i^Hii)'^ is the ionised IGM clumping 
factor. According to the adopted definition of (fesc), Cuu = 1 for a homogeneous 
IGM, or 1 < Ch n ^ 10 taking into account I GM deiisity fl uctuations producing 
the Lya forest (iMiraJda-EscuHe et a].ir200nl : lGneHiTill2nonl 'l. The estimate from 



eq. is rather rough because I have implicitly assumed that the SFR is constant 
during the period of GC formation Atgc- A more realistic SFR as a function of 
redshift requires assuming a specific model for the formation of GCs. I try to 
address this question in the next section. 



4.2. Using the Press- Schechter formalism 

I assume that the formation rate of stars or GCs in galaxies is proportional 
to the merger rate of galaxy halos (each galaxy undergoes a major star burst 
episode when it virializes). Using the Press-Schechter formalism I calculate. 



dt 

dijjl{z) 
dt 



A I dlnMdm, 

Ml dt 



dt 



(3) 
(4) 
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Figure 1. {left) The thin sohd Une in the top panel shows the comoving SFR 
of galaxies as a function of time in our model. The thick solid, dashed and 
short-dashed lines show the SFR of GCs for cases (i) , (ii) and (iii) , respectively 
(assuming fdi = 2). The bottom panel shows the stellar mass budget (in units 
of the baryon abundance), uj( , as a function of time. The segment with arrows 
is a visual aid to compare the GC contribution, cj^^, (assuming fdi = 20) to 

uj( around z ~ 6. (right) Emissivity (photons per baryon per Hubble time) 
as a function of redshift. The thick solid, dashed and short-dashed lines show 
the contribution of GCs for cases (i), (ii) and (iii), respectively (assuming 
fdi = 2). The thin lines show the contribution of galaxies assuming a realistic 
(/esc) = 0.1 X exp[— z/2] (solid) and a constant (fesc) = 5% (dashed). [Plots 
from Ricotti (2002).] 



where Q{Mdm, z)dln{Mdm) is the mass fraction in viriahzed dark matter halos 
of mass Mdm at redshift z. I determine the constants A and by integrating 
eqs. (EI)-® with respect to time, and assuming cogc = 0.1% {i.e., fdi = 2) and 

ioi = l.Auj^: = 13% at z = (the factor 1.4 takes into account the mass loss due 
to stellar winds and SN explosions adopting a continuous star formation law). 
I assume that GCs form in halos with masses Mi < Mdm < The choice 
of Ml and M2 determine the mean redshift, Zf, and time period, Atgc, for the 
formation of old GCs. In order to be consistent with observations I consider 
three cases: case (i) halos with virial temperature 2 x 10^ < Tyir < 5 x 10^ K; 
case (ii) 5 x 10'^ < T^r < 10^ K; and case (iii) 10^ < T^r < 5 x 10^ K. In case 
(i),(ii) and (iii) Atgc = 2.2,3.7 and 5.2 Gyr, respectively, and the GC formation 
rate has a peak at z = 7.5, 6 and 4.6, respectively. Disk and spheroid stars form 
in halos with Mdm > Mm- At 2; > 10 I assume that the first objects form in 
halos with Mm corresponding to a halo virial temperatur e T„ir = 5 x 10^ K. A t 
2; < 10 only objects with T^r > 2 x 10*^ K can form (see lRicotti et al.ll2002lJ ). 

The comoving star formation rate, given by = 'pLol, where p = 5.51 x 10^ Mq 



^By definition il{Mdm., z)dln{Marn) = 1. I find the following values of the constants A = 
1.3%, 1.6%, 0.6% for cases (i), (ii) and (iii) respectively (see text) and B — 12%. 
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Mpc ^ is the mean baryon d ensity at z = 0, i s show n in Fig. ^left). The points 
show the observed SFR from iLanzetta et~ 

In Fig. fright) I show Afph for GCs (thick Hnes) and for galaxies (thin hnes) 
defined as, 



^11 = r^f^^^tn{z), (5) 



duji 



Kh = ri{fesc)^tH{z). (6) 

The thick sohd, dashed and short-dashed hnes show A/^^ for case (i), (ii) and 
(iii), respectively. For comparison, I show (thin solid line) J\f*^ a ssuming (fesc) = 
0.1 X exp[— z/2], derived theoretically bv lRicotti fc Shulll ( 2nnnl'l and normalised 



to fit the observed values (squares) oi J\f*^ at z = 2, 3,4 ( Miralda-Escude et al.l 
I2OOOI ) . The thin dashed line shows N*^ assuming constant {fesc) = 5%. 



5. Conclusions 

The observed Lyman break galaxies at z ~ 3 are probably the most luminous 
starburst galaxies of a population that produced the bulk of the stars in our 
universe. Their formation epoch corresponds to the assembly of the bulges of 
spirals and ellipticals. Nevertheless the observed upper limit on {fesci from 
Lyman break galaxies, (/esc) ~ 10 %, may be insufficient to reionise the IGM 
according to numerical simulations. iFerguson etHI (120021 '). using different ar- 



guments based on the presence of an older stellar population, also noticed that 
the radiation emitted from Lyman break galaxies at z > 3 was insufficient to 
reionise the IGM assuming a continuous star formation mode. 

I propose that GCs during their formation may have produced enough ion- 
ising photons to reionise the IGM. Assuming f^i = 2 [i.e., during their evolution 
GCs have lost half of their original mass), I find a stellar mass fraction in GCs, 
ujgc ~ 0.1%, small compared to the total stellar budget ojI ~ 10% at z = 0. 
But GCs are about 12-13 Gyr old and, if they formed between 5 < 2; < 7 (in 
about 0.5 Gyr), the expected total ojI formed during this time period is about 



ijj 



f 



1%, only 10 times larger than u>gc. Assuming f^i = 20, expected from the 

results of N-body simulations, I find u>gc ~ 1%, suggesting that GC formation is 
an important mode of star formation at high-redshift. The star formation mode 
required to explain the formation of GCs suggests an {fesc) ~ 1 from these ob- 
jects. This is because the mean {fesc) is dominated by the most luminous OB 
associations and GCs are extremely luminous, emitting 5 ~ 3 x 10^^ s~^ ionising 
photons in bursts lasting only 4 Myrs. Moreover, according to many models, 
GCs form in the hot, collisionally-ionised galaxy halo, from which all the ionising 
radiation emitted can escape into the IGM. Therefore it is not too surprising, if 
GCs started forming before z = 6, that their expected contribution to reionisa- 
tion is significant. I find that the number of ionising photons per baryon emitted 
in a Hubble time at 2; = 6 by GCs is A/^^ = {5.1~^^'2)fdi/'^tgc > 1, therefore 
sufficient to reionise the IGM even if we assume fdi = 1. Here, Atgc ~ 0.5 — 2 
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is the period of formation of the bulk of old GCs in Gyrs. Using simple cal- 
culations based on Press-Schechter formalism [see Fig. fright)] I find that, if 
galaxies have {fesc) ^ 5%, GC contribution to reionisation is important. If GCs 
formed by thermal instability in the halo of T^ir ~ 10^ K galaxies (case (iii)), 
the ionising sources have a large bias {i.e., they form in rare peaks of the ini- 
tial density field). Therefore, the mean size of intergalactic Hii regions before 
overlap is large and reionisation is inhomogeneous on large scales. 
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